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Scheme for optimal spin state estimation is considered in 
analogy with phase detection in interferometry. Recently re- 
ported coherent measurements yielding the average fidelity 
(JV+l)/(JV+2) for N particle system corresponds to the stan- 
dard limit of phase resolution 1 / %/N. It provides the bound for 
incoherent measurements when each particle is detected sep- 
arately and information is used optimally. For specific states, 
improvement up to the value 1/N is possible in quantum the- 
ory. The best results are obtained combining sequentially 
coherent measurements on fractional groups of particles. 



I. INTRODUCTION 



measurement on the Stern-Gerlach apparatus and inter- 
ferometry. Particularly, it will be demonstrated that the 
above mentioned score represents the ultimate limitation 
for sequential measurements performed on each particle 
separately for any quantum state. This corresponds to 
the standard resolution 1/yN currently reached in in- 
terferometry. However, this regime does not represent 
the ultimate strategy. In analogy with quantum interfer- 
ometry performance may achieve the resolution up to 
1/N provided that spin orientation is properly coded 
into a quantum state of N particles. In the particular 
case addressed in this contribution the optimal strategy 
corresponds to sequentially performed coherent measure- 
ments. This explicit example demonstrates the complex- 
ity of optimal measurement, which can combine advan- 
tages of both the coherent and sequential measurements 
with groups of particles. 
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Recently Peres and Wotters jjj formulated a conjec- 
ture: Coherent measurement performed on the collective 
system is more efficient than sequential measurements of 
individual particles. This idea has been further developed 
by Massar and Popescu They formulated this con- 
jecture as a proposal for a "quantum game." The player 
has N identical copies of 1/2 spin particles prepared in 
an unknown pure state \ip), and he is allowed to do any 
measurement. Possible results of the measurement will 
be denoted by an index r. The aim of the measurement is 
to determine the original state of the system. Therefore, 
in the next step the measured data should be attached to 
a state |-0) r , which represents the players's estimation of 
an unknown true state. In the last stage of the game the 
true state is compared with its estimate and their coin- 
cidence is quantified by the so-called fidelity: \(ip\tp} r \ 2 . 
The runs are repeated many times with varying true state 
\y!j). The final score of the quantum game is given by the 
averaged fidelity 



s = {\(m r \ 2 }c 



(i) 



where averaging is carried out over the measured data 
and all the true states {} av = {r, Massar and 

Popescu proved that the maximum score is (N+1)/(N+ 
2). This value cannot be reached by measurements act- 
ing on isolated particles. Derka, Buzek and Ekert ||] 
showed that this score can be obtained by coherent mea- 
surement described by a finite-dimensional probability 
operator valued measure. 

The aim of this contribution is to address the rela- 
tion between recently optimized measurement, repeated 



II. ADAPTIVE STERN-GERLACH SPIN 
DETECTION 

Assume a standard measurement on an ideal Stern- 
Gerlach apparatus. A sample particle is prepared in the 
pure spin state 



|n)(n| = i[l + (na)], 



(2) 



where n represents the unity vector on the Poincare 
sphere, (ncr) being its scalar product with vector of Pauli 
matrices. The impinging particles will be deviated up or 
down. In the long run of repetitions the relative fre- 
quencies will approach the prediction of quantum theory. 
Representing the setting of the Stern-Gerlach apparatus 
by the unity vector m, the probabilities of detection "up" 
(+) and "down" (— ) read 



p± = -[l±(mn)]. 



(3) 



What is the best possible but still feasible result, which 
would predict the spin orientation with the highest accu- 
racy? The most accurate state estimation may be done 
if all tested particles were registered in the same output 
channel of the SG apparatus. In such a case the best 
estimation of the spin corresponds to the orientation of 
the SG apparatus. Of course, it does not mean that the 
estimated direction will fit the spin orientation exactly. 
Deviations are distributed according to the posterior dis- 
tribution conditioned by the detected data. This can be 
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handled analogously to the phase estimation |4j . The de- 
viations between the estimated and the true directions 
are given by the Bayes theorem as the posterior proba- 
bility density 

P(n) = ^±±eos 2 >/2) (4) 

over the Poincare sphere. The vector n is parametrized 
by the Euler angles 9, <f> in the coordinate system where 
the direction of z-axis coincide always with the estimated 
direction (i.e. with the orientation of SG apparatus m). 
Notice that this is in accordance with the rules of the 
quantum game as defined in 0]. The result of the mea- 
surement is always a specific direction, namely the setting 
of SG apparatus. In this case the score reads 

S = J P(n)d 2 n n cos 2 (6/2) = ^±1. (5) 

This is the upper bound of sequential measurements on 
single particles with felicitously rotated SG apparatus. 
This is obviously an ultimate limitation since such re- 
sults are possible and none measurement performed se- 
quentially on single particles can yield better spin pre- 
diction. However, on the contrary to the coherent mea- 
surements, this resolution cannot be really achieved, but 
it may be approximated with an arbitrary accuracy for 
N large enough. 

The possible realization may be suggested as an adap- 
tive scheme, where the orientation of SG device depends 
on the previous results. The aim of the scheme is to find 
such an orientation, where almost all the particles are 
counted on the same port. This is obviously always a lit- 
tle worse than ideal case, since some portion of counted 
particles must always be used for corrections of the SG- 
device orientation and are therefore "lost" . Differences 
between ideal and realistic scheme seem to be negligible 
as demonstrated in Fig. [I]. An adaptive measurement 
is simulated here. The procedure starts with project- 
ing a single particle into the three orthogonal randomly 
chosen directions. The choice of the orientation of the 
subsequent SG measurement represents its own interest- 
ing optimization problem. Obviously, for getting the best 
score, it seems to be advantageous to project the particles 
into the most probable orientation. However, this will not 
reveal new corrections to the orientation of SG apparatus 
well. On the other hand, the particles may be counted 
with the SG device oriented perpendicularly to the most 
probable orientation. In this case the spreading of the 
data will be obviously broader than in the former case, 
but such measurement will be more sensitive to the devi- 
ation from the true direction. In the simulation scheme, 
the latter approach has been used and algorithm for syn- 
thesis of incompatible spin projections has been used ||. 
The treatment is still not optimal, however, as will be 
seen in the next paragraph, the differences seem to be 
unimportant. 
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FIG. 1. Numerical simulation of one-particle adaptive mea- 
surement. 



III. ANALOGY BETWEEN SPIN ESTIMATION 
AND INTERFEROMETRY 

There is a clear analogy between spin measurement 
and phase interferometry. As the proper resolution mea- 
sure, the dispersion of phase may be defined as (6|j^] 

D 2 = l-{cos9} 2 av . (6) 

Usually, the averaging is done over the data only. Pro- 
vided that the width of phase distribution is small, dis- 
persion tends to the standard variance of phase. Using 
the definition of the score 

S=l + l{cose} av , (7) 

both the measures fulfil the relation 

S(S + 1) = ^-. (8) 

Consequently, the value S = (N + 1)/(N + 2) is noth- 
ing else as dispersion (phase variance) D«2/ y/N. This 
is the so-called standard limit of phase resolution 
Any standard measurement is scaled in this way and all 
the measurements differ by some multiplicative factors 
only. Loosely speaking, all the classical strategies are 
essentially equivalent from this viewpoint. This is why 
it has perhaps little sense to optimize further the adap- 
tive scheme. For example, provided that one will use the 
most straightforward method of spin estimation based 
on the measurement of x, y, and z components of the 
spin on the Poincare sphere always with N/3 particles, 
the resulting score may be evaluated asymptotically as 
S w 1 - 114/(100iV). The difference between optimal co- 
herent and realistic sequential measurements is assumed 
sometimes to be significant for small number of particles. 
However, in this case all the predictions are rather un- 
certain. In the case of phase detection it has little sense 
to compare two phase distributions, whose widths are 
comparable to 27r window. Then the phase knowledge 
is almost equally bad. Conventional resolution measures 
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possess good meaning only when the effective width is 
substantially less than the width of the interval. 

More profound analogy between the spin measurement 
and interferometry follows from the common nature of 
the SU(2) symmetry H. As well known, the resolution 
up to the order 1/N may be achieved in interferome- 
try. In this case, however, it is recommended to modify 
slightly the proposed quantum game. Suppose that some- 
body wants to communicate an orientation of the axis in 
3D space. Just an axis, not the direction of an "arrow". 
For this purpose N 1/2-spin particles are available and 
any measurement on these particles is allowed. Instead 
of N identical copies, one may consider a general quan- 
tum state spanned by N particles, in which the unknown 
orientation is coded. As the result of the measurement 
the unknown axis should be find. The score is defined 
in the same manner as before. The question are: "How 
to encode information on axis orientation into a quan- 
tum state of the particles? What measurement should 
be done in order to obtain the best score?" 



IV. SUPERRESOLUTION IN 
INTERFEROMETRY AND SPIN ESTIMATION 

Let us review briefly the description of SU(2) interfer- 
ometers. The transformation of an internal state is given 
by an unitary transformation of an input state 



U(0,(p)=exp\-i6J a ((p) = 
exp(z</? J3) exp(— i9 J2) exp(— i<pJz) 



(9) 



where Ji, J2, J3 corresponds to generators of SU(2) group 
The transformation is given by the unity vector 
n = (cos <p sin 9, sin <p sin 9, cos 0) . An input state |in) may 
be any N particle state. The measurement may be rep- 
resented by projectors into an output state |out). The 
posterior phase distribution corresponds to the scatter- 
ing amplitude 



(10) 



This scheme encompasses the Massar, Popescu quan- 
tum game as a special case for the choice of the input 
state I in) = |JV, iV). Score depends on the accuracy of 
detection of 9 variable. As well known, the highest ac- 
curacy is achieved when the phase shift near the zero 
value is detected. This corresponds to the detection of 
the same quantum state on the output |out) = | JV, iV) . 
The particles feeding a single input port of an interfer- 
ometer appears again on the corresponding output port. 
Obviously, the interpretation does not change provided 
that particles will enter and leave the interferometer se- 
quentially. This is the consequence of the famous Dirac 
statement that "each particle interferes with itself." |lC| ] 
Ultimate score (N + l)/(N + 2) is relevant just to this 
regime. Considerations related to adaptive SG detection 
are obviously related also to this situation. 



In the following, for given input state and for given 
setting of the parameters of the interferometer, only the 
most favourable but still feasible output will be assumed. 
Such felicitous results provide an upper bound condi- 
tioned by the realistic measurements. Obviously, no 
other measurement of the given type, can provide bet- 
ter prediction. 

As well known in interferometry |^fl|Tl|| , better resolu- 
tion may be obtained provided that both the input ports 
of an interferometer are fed simultaneously by an equal 
number of r particles. In the case when the same output 
state appears on the output, the phase shift prediction 
will be the sharpest. This corresponds to the input and 
output states |in) = \2r, 0) and |out) = |2r, 0). The scat- 
tering amplitude then depends only on the variable 6 as 

P{9) oc [P°(cos9)] 2 . (11) 

The Legendre function P®(cos9) may be approximated 
by Bessel function Jo(r9) for large index r. Consequently, 
the probability distribution is not integrable for r — > 00 
and must be therefore treated more carefully p2l . Partic- 
ularly, it is therefore not the best strategy to use all the 
energy of N states in the single coherent measurement. 
The particles should be divided into several groups and 
the measurement should be sequentially repeated sev- 
eral times. The accumulation of information is expressed 
mathematically by the multiplication of corresponding 
distribution functions (p]). This tends to narrowing 
of the posterior probability distribution. The optimal 
regime for the 9 detection has been roughly estimated in 
Rcf . [O . Optimal number of repetitions of the coherent 
measurement was expected to be approximately n rs 4. 

Let us interpret this interferometric measurement in 
terms of the spin estimation. The significant difference 
between spin 1/2 particles and photons is connected with 
their fermionic and bosonic nature. The input wave 
function constructed from fermions must be "artificially" 
symmetrized with respect to all the distinguishable par- 
ticles. The input state | JV, 0) is a superposition involving 
all the combinations of N/2 particles with the spin up 
and N/2 particles with the spin down. This N particle 
state is entangled. 
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FIG. 2. The posterior probability distribution provided 
that the exactly same number of spins "up" and "down" has 
been detected. The total number of particles is 20. 
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Let us analyse in detail the the score for such states. 
As before the most favourable but still physically feasi- 
ble situation will be interpreted as an ultimate bound of 
the state detection. Due to the symmetry of the prob- 
lem, the best resolution may be inferred provided that 
the state \N, 0) feeds the SG apparatus. The results of 
such measurement will be deterministic in the case when 
the apparatus is oriented along the spin state of each 
particle. In this case a half of particles will be detected 
with the spin "up" and a half of particles with the spin 
"down" . Assume now that this really happened. What 
can be inferred from this event? Obviously, this might 
occure also in the case when the setting of SG differs from 
the right spin orientation by an angle 8. The probability 
that this happens is proportional to the scattering ampli- 
tude ([n]). The probability is sketched in Fig. ^ - as can 
be seen it shows oscillations. The domain is restricted to 
the values (0,7r/2) since the method just finds the axis 
but not its vector orientation. The score is plotted in 
Fig. [| As shown it does not improve with increasing 
number of particles. It is caused by the heavy tails of the 
posterior distribution (Fig. |^), which does not contain 
the dominant amount of the probability in the central 
peak. Nevertheless, there is a way how to suppress this 
undesired behaviour. Provided that the measurement is 
repeated, the corresponding posterior distribution will be 
given by the product of partial results. Again, the most 
favourable situation for spin estimation is characterized 
by the repetition of the same "optimal" result. This is 
again feasible, provided that the true spin orientation 
and projection do not differ significantly. The resulting 
score in dependence on the total number of particles is 
sketched in Fig. ||. Notice, that SG measurement must 
be done with N/2 particle state in this case, since the 
measurement is repeated twice. This procedure may be 
further generalized including an arbitrary number of rep- 
etitions. As shown in Fig ||, the score increases up to the 
5 repetitions and then starts to decrease. This seems to 
be in good agreement with the rough asymptotical anal- 
ysis presented in Ref. fl^| , where the optimal repetition 
rate has been found as 4. 

This result provides an ultimate bound in the following 
sense: If the given input state is measured with the help 
of SG projections, the score cannot be better than this 
ultimate value. The question whether this value may be 
attained is not answered here. Intuitively, in the case of 
large number of particles it might be well approximated 
by an adaptive scheme as in previous case of standard 
resolution. However, the adaptive scheme cannot be ap- 
plied to the first measurement. The results may therefore 
appear as overestimated in this case. This explains why 
the analysis applied here provides the score S w 0.875 for 
the state | |, J.), whereas the result of Gisin and Popescu 
|L5| gives the value 0.789. In the asymptotical limit, the 
dispersion is approximately given by the relation 



The optimal score reads asymptotically 
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FIG. 3. Ultimate values of score for several repetitions of 
coherent measurements on groups of particles. 

One may ask whether the required states containing 
half particles with spins "up" and half with spins "down" 
(in an arbitrary direction) can be prepared from a set 
of particles with all spins "up" , i.e. if it is possible to 
"turn" arbitrary quantum state into a state orthogonal 
to it. The general answer is no - the linearity of quantum 
mechanics does not allow this. Particularly in case of 
1/2-spin particles, if one is able to do: | |) — ► | J.) and 
| I) — > e llp \ |) then for an arbitrary state it follows from 
linearity of quantum evolution that 

m = (a\ T) + /3| l»-(a| l) + e^| T»< 

The last state is orthogonal to | \P) if and only if arg(a) + 
7r = arg(/3) + (f>. Thus in special situations (correspond- 
ing to "real" subspaces) the mentioned transformation is 
possible (e.g. when spin projections lay in a given plane, 
also for linear polarizations of photon or for interferom- 
etry with fixed splitting ratios and varied a phase differ- 
ence only). But an arbitrary spin cannot be turned to 
the perfectly opposite one. 



V. CONCLUSIONS 

The profound analogy between interferometry and spin 
estimation has been addressed. The performance of both 
the schemes discussed has been conditioned by the re- 
alistic measurements only. As demonstrated, the re- 
cently reported optimal spin estimation corresponds to 
the standard quantum limit characterized by the reso- 
lution 1/y/N. It may be achieved by coherent measure- 
ments and well approximated by sequential ones. Beyond 
this regime, the quantum theory admits the resolution up 
to 1/N. Then, however, quantum interferences must be 
employed. In the case of a thought experiment with an 
ensemble of spin 1/2 particles it requires an entangled 
input state of N particles. Besides this, optimal SG de- 
tection must combine advantages of both the coherent 
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and sequential measurements. This example illustrates 
complexity of the optimal treatment in estimation prob- 
lems. 
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